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Abstract 

We consider the effect of a small cut-off e on the velocity of a traveling wave in one dimension. Simulations done over 
more than ten orders of magnitude as well as a simple theoretical argument indicate that the effect of the cut-off e 
is to select a single velocity which converges when e — > to the one predicted by the marginal stability argument. 
For small e, the shift in velocity has the form K(\oge)~ 2 and our prediction for the constant K agrees very well with 
the results of our simulations. A very similar logarithmic shift appears in more complicated situations, in particular 
in finite size effects of some microscopic stochastic systems. Our theoretical approach can also be extended to give a 
simple way of deriving the shift in position due to initial conditions in the Fisher-Kolmogorov or similar equations. 

PACS: 02.50.Ey, 03.40.Kf, 47.20.Ky 



1 Introduction 

Equations describing the propagation of a front be- 
tween a stable and an unstable state appear jl], g 
^, ^, ^, Q in a large variety of situations in physics, 
chemistry and biology. One of the simplest equations 
of this kind is the Fisher-Kolmogorov jl], ^ equation 



dh 



d 2 h 

dx 2 



+ h-h 3 



(1) 



which describes the evolution of a space and time de- 
pendent concentration h(x,t) in a reaction-diffusion 
system. This equation, originally introduced to study 
the spread of advantageous genes in a population Q , 
has been widely used in other contexts, in particular 
to describe the time dependence of the concentration 
of some species in a chemical reaction ^] . 

For such an equation, the uniform solutions h = 1 
and h = are respectively stable and unstable and it 
is known ]|, ^, O, |l^| that for initial conditions 
such that h(x, 0) — ► 1 as x — > — oo and h[x, 0) — * as 
x — > +oo there exists a one parameter family F v of 
traveling wave solutions (indexed by their velocity v) 



of the form 



h(x, i) = F v (x — vt) , 



(2) 



with F v decreasing, F v (z) — > 1 as z — > —oo and 
F v (z) — > as z — > oo. The analytic expression of 
the shape F v is in general not known but one can de- 
termine the range of velocities v for which solutions 
of type (||) exist. If one assumes an exponential decay 



F v (z) 



for large z, 



(3) 



it is easy to see by replacing 
the velocity v is given by 



and (g) into (Q) that 



u(7) =7 + -. 

7 



(4) 



As 7 is arbitrary, this shows the well known fact 
that the range of possible velocities is v > 2. The 
minimal velocity vq = 2 is reached for 70 — 1 and for 
steep enough initial conditions h(x, 0) (which decay 
faster than e~^° x ), the solution selected [|[ § || [f| 



_ L0[ , |11| , |12J for large t is the one corresponding to this 
minimal velocity vq. 
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Equations of type (11) are obtained either as the 
large scale limit |5[ |[ 13|, [l4|, [l5| |l6|] or as the mean 
field limit jl?]] of physical situations which are dis- 
crete at the microscopic level (particles, lattice mod- 
els, etc.) As the number of particles is an integer, the 
concentration h(x, t) could be thought as being larger 
than some e, which would correspond to the value of 
h(x, t) when a single particle is present. Equations of 
type (|l|) appear then as the limit of the discrete model 
when e — > 0. Several authors ||, [l3|, [l4| have already 
noticed in their numerical works that the speed v e of 
the discrete model converges slowly, as e tends to 0, 
towards the minimal velocity vq. We believe that the 
main effect of having e 7^ is to introduce a cut-off in 
the tail of the front, and that this changes noticeably 
the speed. 

The speed of the front is in general governed by 
its tail. In the present work, we consider equations 
similar to (|l|), which we modify in such a way that 
whenever h(x, t) is much smaller than a cut-off e, it 
is replaced by 0. The cut-off e can be introduced by 
replacing (Q) by 



with 



a{h) 
a(h) 



= 1 
< 1 



if h > e, 
if h < e. 



(5) 



(6) 



For example, one could choose a(h) = 1 for h > e and 
a(h) = h/e for h < e. Another choice that we will use 
in section ^ is simply a(h) = 1 if h > £ and a(h) = 
if h < e. 

The question we address here is the effect of the 
cut-off e on the velocity v e of the front. We will 
show that the velocity v £ converges, as e — ► 0, to the 
minimal velocity vq of the original problem (without 
cut-off) and that the main correction to the velocity 
of the front is 



v"(lo) 



1 



(loge)' 



(7) 



for an equation of type (|TJ) for which the velocity is 
related to the exponential decay 7 of the shape W) 



by some relation ^(7). (Everywhere we note by vo 
the minimal velocity and 70 the corresponding value 
of the decay 7.) In the particular case of equation 
([!]), where ^(7) is given by (^), this becomes 



(8) 



(loge)' 



In section || we describe an equation of type (Q) 
where both space and time are discrete, so that simu- 
lations are much easier to perform. The results of the 
numerical simulations of this equation are described 
in section ||: as e — > 0, the velocity is seen to con- 
verge like (logs) 2 to the minimal velocity i>o, and 
the shape of the front appears to take a scaling form. 

In section || we show that for equations of type ([!]) 
in presence of a small cut-off e as in (||), one can 
calculate both the shape of the front and the shift in 
velocity. The results are in excellent agreement with 
the numerical data of section 0. 

In section^ we consider a model defined, for a finite 
number N of particles, by some microscopic stochas- 
tic dynamics which reduces to the front equation of 
sections || and [I] in the limit N — > 00. Despite the 
presence of noise, our simulations indicate that in this 
case too, the velocity dependence of the front decays 



slowly (as (logiV)' 
the front. 



to the minimal velocity vq of 



2 A discrete front equation 

To perform numerical simulations, it is much easier to 
study a case where both time and space are discrete 
variables. We consider here the equation 

h{x,t+l)=g(x,t) 0[g(x,t)-e], (9a) 

where 

g(x, t) = 1 - [1 - ph{x - 1, t) - (1 - p)h(x, t)] 2 . 

(9b) 

Time is a discrete variable and if initially the concen- 
tration h(x, 0) is only defined when x is an integer, 
h(x,t) remains so at any later time. Because t and x 
are both integers, the cut-off e can be introduced as 



2 



in (0) in the crudest way using a Heaviside Q func- 
tion. (Wc have checked however that other ways of 
introducing the cut-off e as in <M, ^) do not change 
the results.) 

Equation (||) appears naturally (in the limit e = 0) 
in the problem of directed polymers on disordered 
trees @ (where the energy of the bonds is ei- 
ther 1 with probability p or with probability 1 — p) . 
At this stage we will not give a justification for in- 
troducing the cut-off £. This will be discussed in 
section ||. 

We consider for the initial condition a step function 

h(x,0) =0 ifx>0, (10) 

h(x,0) =1 ifx<0. 

Clearly for such an initial condition, h(x, t) = 1 for 
x < at all times. As h(x, t) ~ 1 behind the front 
and h(x,t) ~ ahead of the front, we define the 
position X t of the front at time t by 

^ = 5^/1(1,*). (ii) 

The velocity of the front v e can then be calculated 

by 

v, = lim ^ = (AVi - A t ) , (12) 

where the average is taken over time. (Note that 
as h(x,t) is only defined on integers, the difference 
Xt+i — X t is time dependent and has to be averaged 
as in @.) 

When e — 0, the evolution equation (0) becomes 

h(x, t + 1) = 1 - [l - ph(x — 1, t) — (1 — p)/i(ic, t)] 2 . 

(13) 

As for (|lj) , there is a one parameter family of solutions 
F v of the form (||) indexed by the velocity v which is 
related (||) to the exponential decay 7 of the shape 

by 

«( 7 ) = i log(2pe 7 + 2(1 -p)). (14) 



(This relation is obtained as (Q) by considering the 
tail of the front where h(x,t) is small and where there- 
fore @ can be linearized.) 

One can show that for p < i, ^(7) reaches a mini- 
mal value Vq smaller than 1 for some 70, whereas for 
P > hi v (l) is a strictly decreasing function of 7, im- 
plying that the minimal velocity is vq = lim ^(7) = 

7 — >oc 

1. 

We will not discuss here this phase transition and 
we assume from now on that p < |. Table [l] gives 
some values of vq and 70 obtained from (|l4|). 



p 


0.05 


0.25 


0.45 


7o 


2.751111 ... 


2.553 244 ... 


4.051851 ... 


vo 


0.451818 ... 


0.810710 ... 


0.979187 ... 



Table 1: Values of 70 and vq for some p when e = 0. 



It is important to notice that for p < i , the func- 
tion 11(7) has a single minimum at 70. Therefore, 
there are in general two choices 71 and 72 of 7 for 
each velocity v. For v ^ the exponential decay of 
F v (z) is dominated by min(7i,72)- As v — > vq, the 
two roots 71 and 72 become equal, and the effect of 
this degeneracy gives (in a well chosen frame) 

F V0 (z)~Az e~ 70Z for large z, (15) 

where A is a constant. This large z behavior can be 
recovered by looking at the general solution of the 
linearized form of equation ( |l3| ) 

h(x, t + l) = 2ph(x - 1, t) + 2(1 - p)h(x, t). (16) 

3 Numerical determination of 
the velocity 

We iterated numerically (J9J) with the initial condition 
( |To| ) for several choices of p < \ and for e varying 
between 0.03 and 10 -17 . We observed that the speed 
is usually very easy to measure because, after a short 
transient time, the system reaches a periodic regime 
for which 

h{x, t + T) = h(x — Y, t) (17) 



3 



for some constants T and Y . 
front is then simply given by 

Y 



The speed v £ of the a given time t relative to its position X t by 



(18) 



For example, for p — 0.25 and e = 10 5 , we find 
T = 431 and Y = 343 so that v e = 343/431. The 
emergence of this periodic behavior is due to the lock- 
ing of the dynamical system of the h(x,t) on a limit 
cycle. Because Y and T are integers, our numerical 
simulations give the speed with an infinite accuracy. 

For each choice of p and e. we measured the speed 
of the front, as defined by (|l2| ) and its shape. Figure [j] 
is a log-log plot of the difference v — v e versus e 
(varying between 0.03 to 10~ 17 ) for three choices of 
the parameter p. The solid lines on the plot indicate 
the value predicted by the calculations of section |[ 
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Figure 1: The difference vq — v £ for p = 0.05, 0.25 
and 0.5. The symbols represent the result of our nu- 
merical simulations and the solid lines indicate the 
prediction of the analysis of section [|. 

We see on this figure that the velocity v E converges 
slowly towards the minimal velocity Vq as e — > 0. 
Our simulations, done over several orders of magni- 
tude (here, fifteen), reveal that the convergence is 
logarithmic: vq — v e ~ (loge) 

As the front is moving, to measure its shape, we 
need to locate its position. Here we use expression 
(O) and we measure the shape s e {z) of the front at 



,(z) = h(z + X t ,t) 



(19) 



When the system reaches the limit cycle (p7|), the 
shape s £ (z) becomes roughly independent of the time 
chosen. (In fact it becomes periodic of period T, but 
the shape s e has a smooth envelope.) We have mea- 
sured this shape at some arbitrary large enough time 
to avoid transient effects. As we expect s e (z) to look 
more and more like F VQ (z) as e tends to 0, we nor- 
malize this shape by dividing it by e~ loZ . The re- 
sult s e (z)e loZ is plotted versus z for p = 0.25 and 
e = 10~ 9 , 10" 11 , 10~ 13 , 10" 15 and 10" 17 in figure |. 




10 12 14 16 



Figure 2: Normalized shape of the front s e (z)e y ° z 
versus z for p = 0.25 and several choices of e. 

On the left part of the graph, our data coincide 
over an increasing range as e decreases, indicating 
that far from the cut-off, the shape converges to ex- 
pression (|l^) of F Vq (z). On the right part, the curves 
increase up to a maximum before falling down to 
some small value which seems to be independent of s. 
When e is multiplied by a constant factor (here 10~ 2 ), 
the maximum as well as the right part of the curves 
are translated by a constant amount. This indicates 
that for e small enough, the shape s £ (z) in the tail 
(that is for z large) takes the scaling form 



s e (z)~\\oge\ G 



|loge| 



(20) 
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We will see that our analysis of section ^ does pre- 
dict this scaling form. As one expects this shape to 
coincide with the asymptotic form ( |l5| ) of F Vo (z) for 
1 <C z <C | logs], the scaling function G(y) should be 
linear for small y. 

4 Calculation of the velocity for 
a small cut-off 

The first remark we make is that as soon as we in- 
troduce a cut-off through a function a{h) which is 
everywhere smaller than 1, the velocity v £ of the 
front is lowered compared to the velocity obtained 
in the absence of a cut-off. This is easy to check 
by comparing a solution h £ (x,t) of (5) where a(h) 
is present and a solution /z-o (a:, t) of (1). If initially 
h £ (x,0) < ho(x,0), the solution h £ will never be able 
to take over the solution h . Indeed, would the two 
functions h £ (x,0) and ho(x,0) coincide for the first 
time at some point x, we would have at that point 
d 2 h £ /dx 2 < d 2 ho/dx 2 and together with the effect of 
a(h) this would bring back the system in the situa- 
tion where h £ (x,t) < ho(x,t) ||, This shows that 
v £ < v - 

For the calculation of the velocity v £ , we will 
consider first the modified Fisher-Kolmogorov equa- 
tion (|^) when the cut-off function a{h) is simply given 

by 



a{h) = ®{h-e). 



(21) 



In this section we will calculate the leading correc- 
tion to the velocity when e is small and we will obtain 
the scaling function G which appears in (^). Then 
we will discuss briefly how our analysis could be ex- 
tended to more general forms of the cut-off function 
a(h) or to other traveling wave equations such as (Q). 

As v e is the velocity of the front, its shape s £ (z) = 
h{z + v £ t, t) in the asymptotic regime satisfies 

v s s' e + s" + (s e - s e )a(s e ) = 0. 

When e is small, with the choice ( |2l| ) for a(h), we can 
decompose the range of values of z into three regions: 

Region I where s e (z) is not small com- 
pared to 1. 



Region II where e < s e (z) -C 1. 
Region III where s e (z) < e. 

In region I, the shape of the front s £ looks like F Vo 
whereas in regions II and III, as s £ is small, it satisfies 
the linear equations 



v e s' + s'! 



= 



v £ s' + s" = 



in region II, (22) 
in region III. (23) 



These linear equations (^,^3|) can be solved easily. 
The only problem is to make sure that the solution 
in region II and its derivative coincides with _F!„ at 
the boundary between I and II and with the solution 
valid in region III at the boundary between II and III. 
If we call A the shift in the velocity 



A = v - v E 



(24) 



and if we note 7 r ± the two roots of the equation 
v {l) — v ei the shape s £ is given in the three regions 

by 

s £ (z) ~ F Vo (z) in region I, 

s e (z) ~ Ce" 7 ' z sin(7iZ + D) in region II, (25) 



s e (z) 



se 



-v e {z—zo) 



in region III, 



and we can determine the unknown quantities C, D, 
zq and v e by using the boundary conditions. 

For large z we know from ([[5]) that F VQ {z) ~ 
Aze~ lQZ for some A. Therefore, as 70 — 7 r ~ A 
and 7i ~ A 1 / 2 , the boundary conditions between re- 
gions I and II impose, to leading order in A 1 / 2 , that 
C = A/ji and D = 0. 

At the boundary between regions II and III, we 
have s £ (z) — e and z — z$. If we impose the conti- 
nuity of s £ and of its first derivative at this point, we 
get 



Ae lrZ " sin(7iz ) = £7^, 



(26a) 



and 



Ae lrZo [-j r sin(7iZ ) + % cos(7iZ )] = -v £ e^i. 

(26b) 

Taking the ratio between these two relations leads to 



tan(7i2;o) 



(27) 
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When A is small, j r ~ 70 = 1, v e ~ vq = 2 and 
7i ~ A 1 / 2 . Thus the only way to satisfy ( |27| ) is to 
set 7^0 ~ it and 7r — 7^0 — 7i ~ A 1 / 2 . Therefore, 
( p6| ) implies to leading order that z ~ — (loge)/7o 
and the condition 7^0 ~ 7r gives 



7T 7T7o 

7i - — - Tj 1 

z loge 



(28) 



Then, as 7, is small, the difference A = «o — w £ is 
given by 



v 



v e ~ ~w"(7ohi 



^(7o)^ 2 7o 2 
2(loge) 2 



(29) 



which is the result announced in (0) and (g). 

A different cut-off function a(/i) should not affect 
the shape of s £ in the region 11 or the size Zq of re- 
gion II. Only the precise matching between regions II 
and III might be modified and we do not think that 
this would change the leading dependency of zq in 
e which controls everything. In fact there are other 
choices of the cut-off function a(h) (piecewise con- 
stant) for which we could find the explicit solution 
in region III, confirming that the precise form of a(h) 
does not change (|2^) . The generalization of the above 
argument to equations other than (Q) (and in par- 
ticular to the case studied in sections || and ||) is 
straightforward. Only the form of the linear equation 
is changed and the only effect on the final result (0) 
is that one has to use a different function ^(7). 

When expression (0) is compared in figure |l| to the 
results of the simulations, the agreement is excellent. 
Moreover, in region II, one sees from (^) and ( p8f ) 
that 



s £ (z) ~ loge sin T , | e 



(30) 



which also agrees with the scaling form (20). 

Recently, for a simple model of evolution |T^, pof 
governed by a linear equation, the velocity was found 
to be the logarithm of the cut-off to the power |. 
This result was obtained by an analysis which has 
some similarities to the one presented in this section. 



5 A stochastic model 

Many models described by traveling wave equations 
originate from a large scale limit of microscopic 
stochastic models involving a finite number N of par- 
ticles [^|[ [l4| Here we study such a mi- 
croscopic model, the limit of which reduces to fll3| ) 
when N — > 00. Our numerical results, presented be- 
low, indicate a large N correction to the velocity of 
the form vn — f — a (log N) 2 with a coefficient a 
consistent with the one calculated in section || for 

£ = JV' 

The model we consider in this section appears in 
the study of directed polymers [H and is, up to mi- 
nor changes, equivalent to a model describing the dy- 
namics of hard spheres jl5) . It is a stochastic process 
discrete both in time and space with two parameters: 
N, the number of particles, and p, a real number be- 
tween and I. At time t (t is an integer), we have 
N particles on a line at integer positions x±(t), X2(t), 
... , Xisr(t). Several particles may occupy the same 
site. At each time-step, the TV positions evolve in the 
following way: for each i, we choose two particles ji 
and j[ at random among the N particles. (These 
two particles do not need to be different.) Then we 
update Xi(t) by 

Xi(t + 1) = max(x Ji (i) + c^, Xf.(t) + a'A, (31) 

where cti and a[ are two independent random num- 
bers taking the value 1 with probability p or with 
probability 1 — p. The numbers cti, o^, ji and j[ 
change at each time-step. Initially (t = 0), all parti- 
cles are at the origin so that we have Xi(0) — for 
all i. 

At time t, the distribution of the Xi(t) on the line 
can be represented by a function h{x, t) which counts 
the fraction of particles strictly at the right of x. 



h(x, t) 



- y 



1. 



(32) 



Xi (t)>x 



Obviously h(x,t) is always an integral multiple of 
At t = 0, we have h(x, 0) = 1 if x < and h(x, 0) = 
if rr > 0. One can notice that the definition of the 
position Xt of the front used in (Thy) coincides with 
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the average position of the N particles 

+00 N 
h.(x. A = 

N 



+00 ^ N 



(33) 



Given the positions Xi(t) of all the particles (or, 
equivalently given the function h(x,t)), the Xi(t + 
1) become independent random variables. Therefore, 
given h(x, t), the probability for each particle to have 
at time t + 1 a position strictly larger than x is given 

by 

(h(x,t + l)\h(x,t)) (34) 
= 1 - [l - ph(x - 1, t) - (1 - p)h(x, t)] 2 . 

The difficulty of the problem comes from the fact 
that one can only average h(x,t + 1) over a single 
time-step. On the right hand side of ( p4| ) we see 
terms like h 2 (x,t) or h(x — l,t)h(x,t) and one has 
to calculate all the correlations of the h(x, t) in or- 
der to find (h(x,t + 1)). This makes the problem 
very difficult for finite N. However, given h(x,t), the 
Xi(t + 1) are independent and in the limit N — > 00, 
the fluctuations of h(x, t + 1) are negligible. There- 
fore, when N — > 00, h(x, t) evolves according to the 
deterministic equation (|13|). As the initial condition 
is a step function, we expect the front to move, in the 
limit N — > 00, with the minimal velocity v n of (fl4|). 

For large but finite N, we expect the correction to 
the velocity to have two main origins. First, h(x,t) 
takes only values which are integral multiples of 
so that jf plays a role similar to the cut-off e of sec- 
tion ^. Second, h(x, t) fluctuates around its average 
and this has the effect of adding noise to the evo- 
lution equation (|l3|). In the rest of this section we 
present the results of simulations done for large but 
finite iV and we will see that the shift in the velocity 
seems to be very close to the expression of section [| 
when e= jr. 

With the most direct way of simulating the model 
for N finite, it is difficult to study systems of size 
much larger than 10 6 . Here we use a more sophis- 
ticated method allowing N to become huge. Our 
method, which handles many particles at the same 
time, consists in iterating directly h(x,t). 



Knowing the function h(x,t) at time t, we want 
to calculate h(x,t + 1). We call respectively x m i n 
and i m „ the positions of the leftmost and rightmost 
particles at time t and / = x max — x m j n + 1. In terms 
of the function h(x,t), one has < h(x,t) < 1 if 
and only if x mnl < x < x max . Obviously, all the 
positions Xiit + 1) will lie between x m i n and x max + 
1. The probability pk that a given particle i will be 
located at position x m ; n + k at time t + 1 is 

p k = (h(x min + k - l,t+ 1)) - (h(x mhl +k,t+l)), 

(35) 

with (h{x,t + 1)) given by (|34|). Obviously, pk 7^ 
only for < k < I. 

The probability to have, for every k, nk particles 
at location x m i n + k at time t + 1 is given by 



P(rao, nx, . . . , ni) 



n !ni! 



ni 



rin ni ni 

Po Pi ---Pi 



x 5(N — no — n% — ■ ■ ■ — n{). 

(36) 

Using a random number generator for a binomial 
distribution, expression (^6|) allows to generate ran- 
dom rife. This is done by calculating n according to 
the distribution 



P(n ) = 
then ni with 
P(ni I no) 



n \(N-no)\ 

(JV-no)! 



Po°(l-Po) 



N-n Q 



(37) 



Pi 



nil (N - n - nx)\ \l-p 0/ 

N—hq —ni 

" 1 1 ^— I , (38) 



1 - po 

and so on. This method can be iterated to produce 
the l+l numbers no, n\, . . . , n; distributed according 
to (j36"l). Then we construct h(x,t + 1) by 



h(x,t + l) =1 



if X < X m i n , 



h(x,t + l) 



1 

'N 



E 



if. 



< X < Xr, 



and x — ^min k^ 



i=k+l 

h{x, t + 1) =0 if x > x max + 1. 



(39) 
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As the width I of the front is roughly of order log N, 
this method allows N to be very large. 

Using this method with the generator of random 
binomial numbers given in plj , we have measured the 
velocity vn of the front for several choices of p (0.05, 
0.25 and 0.45) and for N ranging from 100 to 10 16 . 
We measured the velocities with the expression 



VN 



X w e — X 



10 5 



910 5 



(40) 



I'N 
J_ 

N ■ 



Figure 51 is a log-log plot of the difference vq 
versus jj compared to the prediction ([?]) for e 
The variation of vn when using longer times or dif- 
ferent random numbers were not larger than the size 
of the symbols. 




10" 



10" 



10" 



10" 



10" 



10" 



10" 



0.01 



Figure 3: The difference vq — vn versus for three 
choices of p. The symbols represent the result of our 
numerical simulations of the stochastic process and 
the solid lines indicate the prediction (0) for e 



j_ 

N ■ 



We see on figure || that the speed vn of the front 
seems to be given for large N by 



K 



V N ~ V 



(log NY 



(41) 



where the coefficient K is not too different from the 
prediction ([?]). 

The agreement is however not perfect. The 
shift vq — vn seems to be proportional to (logiV) 2 , 
but the constant looks on figure [| slightly different 



from the one predicted by (Q). A possible reason for 
this difference could have been the discretization of 
the front: instead of only cutting off the tail as in 
sections || and ||, here the whole front h(x,t) is con- 
strained to take values multiple of -k. One might 
think that this could explain this discrepancy. How- 
ever, we have checked numerically (the results are 
not presented in this paper) that equation (|l3|) with 
h(x, t) constrained to be a multiple of a cut-off e does 
not give results significantly different from the sim- 
pler model of sections || and [| with only a single 
cut-off. So we think that the full discretization of 
the front can not be responsible for a different con- 
stant K. The discrepancy observed in figure || is more 
likely due to the effect of the randomness of the pro- 
cess. It is however not clear whether this mismatch 
would decrease for even larger TV. It would be inter- 
esting to push further the numerical simulations and 
check the TV-dependence of the front velocity for very 
large N. 



6 Conclusion 

We have shown in the present work that a small cut- 
off e in the tail of solutions of traveling wave equa- 
tions has the effect of selecting a single velocity v e 
for the front. This velocity v e converges to the min- 
imal velocity vq when e — ► and the shift vq — v £ is 
surprisingly large (0, ||). 

Very slow convergences to the minimal velocity 
have been observed in a number of cases 
as well as the example of section [| As the effect of 
the cut-off on the velocity is large, it is reasonable 
to think that it would not be much affected by the 
presence of noise. The example of section || shows 
that the cut-off alone gives at least the right order 
of magnitude for the shift and it would certainly be 
interesting to push further the simulations for this 
particular model to see whether the analysis of sec- 
tion |^ should be modified by the noise. The numer- 
ical method used in section to study a very large 
(N ~ 10 16 ) system was very helpful to observe a log- 
arithmic behavior. We did not succeed to check in 
earlier works |l4], [l5| whether the correction 
was logarithmic, mostly because the published data 



8 



were usually too noisy or obtained on a too small 
range of the parameters. Still even if the cut-off was 
giving the main contribution to the shift of the ve- 
locity, other properties would remain very specific to 
the presence of noise like the diffusion of the position 
of the front |§. 

Our approach of section ^ shows that the effect of a 
small cut-off is the existence of a scaling form [30|) 
which describes the change in the shape of the front 
in its steady state. The effect of initial conditions 
for usual traveling wave equations (with no cut-off) 
leads to a very similar scaling form for the change in 
the shape of the front in the transient regime. This 
is explained in the appendix where we show how the 
logarithmic shift of the position of a front due to ini- 
tial conditions fllQ, PJJ] can be recovered. 



A Effect of initial conditions 
on the position and on the 
shape of the front 

In this appendix we show that ideas very similar to 
those developed in section |J allow one to calculate the 
position and the shape at time t of a front evolving 
according to ([]]), or to a similar equation, given its 
initial shape. The main idea is that in the long time 
limit, there is a region of size \ft ahead of the front 
which keeps the memory of the initial condition. We 
will recover in particular the logarithmic shift in the 
position of the front due to the initial condition |]Io| , 
G3| , namely that if the initial shape is a step function 



we will show that for v > — 2 

X t ~2t+^—±\ogt, 



(A.4) 



whereas the shift is given by (A. 2) for v < —2. Here, 
there is no cut-off but the transient behavior in the 
long time limit gives rise to a scaling function very 
similar to the one discussed in section ^. 

If we write the position of the front at time t as 



X t = v t - c(t), 



(A.5) 



we observed numerically (as in figure ^| of section [jj) 
and we are going to see in the following that the shape 
of the front takes for large t the scaling form 



/,(.;-./) = ra ( ^— ^ ) < "" A ' 



(A.6) 



very similar to (|20| , |30| ) . 

If we use (|A.5| ) and ( A.6 ) into the linearized form 
of equation (|lj), we get using, the fact that vq = 2 
and 70 = 1, 

^G" + -l-(azG' -aG)+t a bG = cG', (A.7) 

where z — (x — X t )t~ a . By wri ting that the leading 
orders of the different terms of (A..7) are comparable, 
we sec that we must have 



f 

2' 

£ 
t • 



(A.8) 
(A.9) 



h(x, 0) = if x > 0, (A.l) f or some ^ anc j tnat tne right hand side of < $Jl\) is 

h(x,0) = 1 if x < 0, negligible. Therefore, the equation satisfied by G is 



then the position X t of the front at time t increases 
like 



dzi G+ 2dz G 



G = 0, 



(A.10) 



X t ~2t 



More generally, if initially 



log*. 



h(x,0) 
h(x,0) 



if x > 0, 
if x < 0, 



(A 2) anc ^ ^ e P os ^i° n °f the front is given by 

X t ~ vot-piogt. (A.ll) 
As in section Q we expect that as t — > oo, the front 

(A.3) 

will approach its limiting form and therefore that for 



z small, the shape will look like (15). Therefore we 
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choose the solution Gp{z) of (A. 10) which is linear 
at z — 0. This solution can be written as an infinite 



Gp(z)=A 



{-iy 



n-l 



^ (2n + l)\ Z 

n=0 y ' i=0 



(-1)™ T(n + j3) 



,2n+l 



(2n + l)! 



(A.12) 



(The second expression is not valid when (3 is a non- 
positive integer.) 

To determine (3, one can notice that the scaling 
form ( |A.6| ) has to match the initial condition when x 
is large and t of order 1. We thus need to calculate 
the asymptotic behavior of G(z) when z is large. 

For certain values of (3, there exist closed expres- 



sions of the sum (A.12). For instance 



G- 2 (z 

Gr(z 

2 v 

G-i(* 
G=(z 

2 

G (z 
G 3 (z 

Gx(z 
Gi(z 



z 

60 
£ 
GO 



A (* + T 

A( z -f! 

a, ,4 



A 



e 4 



(A.13) 



Ae ~ 
A 



e 4 dt, 



o 

e 4 at, 



One can check directly on (A.1C) that G/3 has a sym- 
metry 

G p («) = -ie~ z2/i Gi_ {iz). (A.14) 

For any /?, one can obtain the large z behavior of 
G(z). To do so, we note that for f3 > 0, one can 
rewrite ( A.12 ) as 

A t°° - 



r(/3) Jo 

2A 



1-2/3 



eft i 2/3 ~ 2 sin(i)e7J05) 



For < j3 < 1, the second integral in (A. 15) has a 
non zero limit and this gives the asymptotic behavior 
of Gp(z) 

2 A 

Gp{z) ~ cos(7r/3) r(2/3 - ^z 1 ' 20 . (A.16) 



r(/3) 



From ( A.12 ), one can also show that 

r(/? + i), 



(-in 

u f3 



-G 



13+1, 



(A.17) 



implying that ( A.16| ) remains valid for all /3 except 



for (3 = |, |, |, etc., where the amplitude in (A.16) 
vanishes. For these values of /?, Gp(z) decreases faster 
than a power law (see ( A.13| )). 

The functions Gp calculated so far are acceptable 
scaling functions for the shape of the front only for 



j3 < §. Indeed, one can see in ( |A.16| ) that for | < j3 < 



| the function Gp{z) is negative for large z. In fact, 
for all j3 > |, this function changes its sign at least 
once, so that the scaling form ( |A.6| ) is not reachable 
for an initial h(x, 0) which is always positive. It is 
only for (3 < § that Gp remains po sitive for all z > 0. 

Looking at the asymptotic form ( A.16 ), we see that 
if initially h(x, 0) = x v e mx , the only function Gp(z) 
which has the right large z behavior is su ch tha t 1 — 
2(3 = v, a nd t his gives, together with ( A. 11 ), the 
expression ( |A.4[ ) for the shift of the position. As the 
cases j3 > § are not reachable, all initial conditions 
corresponding to v < —2 or steeper (such as step 
functions) give rise to Gz and the shift in position 
given by JA^ ). 

All the analysis of this appendix can be extended 
to other traveling wave equations such as ([l3|), with 
more general functions v (7) (having a non-degenerate 
minimum at 70) as in (|14|). Then the expres- 



sions ( 



and 



of the shift become 



v t - — log* 

27o 



X t ~ v t 



\-v 



log*. 



(A.18) 



(A.19) 



We thank C. Appert, V. Hakim and J.L. Lebowitz 
for useful discussions. 
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